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TWO-DIMENSIONAL HEAT TRANSFER IN RADIATING STAINLESS- STEEL-CLAD COPPER FINS 
by Norbert 0. Stockman and Edward C. Bittner 
Lewis Research Center 


SUMMARY 

An analysis of the two-dimensional heat transfer in a radiating composite 
fin is presented together with details of. the numerical method of solution. 
Several stainless-steel-clad copper fins typical of finned-tube space radiators 
are analyzed over a range (0 to l) of stainless-steel to copper thickness ratio. 
Items discussed are temperature profiles along the fin length, temperature drop 
across fin thickness, effective conductivity, fin efficiency, heat rejection 
rate per unit weight, and ratio of effective conductivity to effective density. 
It was found for the particular geometries analyzed that the two-dimensional 
temperature distribution in the fin cross section is essentially one- 
dimensional. The use of the thickness -averaged conductivity in one-dimensional 
calculations yields essentially the same results as the two-dimensional calcu- 
lations . 


INTRODUCTION 

The radiator of an electric -power generation system for space applications 
may account for a substantial portion of the powerplant total weight (refs. 1 
and 2). An important factor in the minimum-weight design of these radiators is 
the selection of suitable materials (ref. 3). The large number of basic re- 
quirements (see ref. 3 for discussion) that must be satisfied makes the selec- 
tion difficult. Often all the desirable qualities cannot be found in any one 
currently available material, and composite materials are therefore considered. 
In particular, the fin material for a space powerplant radiator should have a 
high value for the ratio of thermal conductivity to density (ref. 3), adequate 
strength, and a low sublimation rate in vacuum at high temperature. Stainless 
steel has adequate strength and satisfactory sublimation rate but low thermal 
conductivity, whereas copper has very high conductivity but questionable 
strength and sublimation rate at temperatures above about 800° K. These facts 
lead to a consideration of the composite material, stainless- steel-clad copper, 
for radiator fins (e.g., in the SNAP-50/SPUR design studies, refs. 4 and 5). 
Stainless-steel-clad copper has already been used in other applications (e.g., 
cooking utensils and heat exchangers) for several years and is commercially 
available. If a composite material is to be used for radiator fins, its effec- 
tive thermal conductivity must be known for use in one -dimensional calculations 
for radiator design and comparison purposes. Suppliers of stainless -steel-clad 
copper material usually quote two values of effective conductivity: one value 
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for one -dimensional conduction in the direction parallel to the plane of the 
interface, and another value for one -dimensional conduction normal to the plane 
of the interface. For any significant thickness of clad these two values are 
quite different (e.g., for a stainless-steel to copper thickness ratio of only 
0.05 these values are, respectively, 0.955 and 0.525 that of copper alone). In 
a radiator fin the conduction heat transfer is two-dimensional, and the value 
of effective conductivity will lie somewhere between the two one-dimensional 
values . 

Specifically, the heat radiated by a fin is initially supplied across the 
fin thickness at the tube edge and is conducted parallel to the fin surface. 

The heat must subsequently reach the surface of the fin where it is radiated to 
space. There is clearly a two-dimensional effect involved, and the magnitude 
of this effect will increase as the stainless-steel to copper thickness ratio 
increases since the stainless steel has a much lower conductivity than copper. 
On the other hand, a radiator fin is usually very long relative to its thick- 
ness, which will tend to decrease the two-dimensional effect. Therefore, it is 
not a priori evident how great the two-dimensional effect will be, nor is it 
evident what the effective conductivity of a radiating fin will be. 

Reference 6 attacks this type of problem in a quasi- two-dimensional way by 
assuming a one -dimensional heat flow in the fin core parallel to the radiating 
surface and one-dimensional heat flow in the clad or coating normal to the ra- 
diating surface. The primary application of the analysis of reference 6 is for 
coated fins having a relatively thin, very low conductivity, nonmetallic coat- 
ing. The metallic-clad and the nonmetallic-coated fins present essentially the 
same problem. However, it is not obvious that the method of reference 6 will 
be adequate for metallic-clad fins in all cases because a metallic clad is 
thicker and has a higher value of thermal conductivity than a nonmetallic coat- 
ing. 


This report presents a brief analysis of the two-dimensional heat transfer 
in a composite fin with the conventional fin-tube radiator boundary conditions, 
that is, radiation from both surfaces and prescribed constant temperature at 
the fin base. An exact effective conductivity is obtained from the two- 
dimensional results and is compared with the conductivity for one -dimensional 
conduction in the direction parallel to the fin surface. 

Several cases are analyzed that cover a range of the conductance parameter 
(the conventional identifying parameter for radiating fins). The range in this 
parameter is obtained by varying the length and thickness of the fin at two 
different temperature levels. For each case several values of stainless-steel 
to copper thickness ratio are prescribed. The results to be presented are 
discussed under two main headings: first, the degree of two- dimensionality of 

the several cases, and second, the effect of varying the conductance parameter 
and the stainless -steel to copper thickness ratio on fin temperature profiles, 
fin efficiency, heat rejected per unit weight, and effective conductivity to 
effective density ratio. 
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(a) Application. 


ANALYSIS 

The configuration to be analyzed 
is shown in figure 1. This is a 
stainless-steel-clad copper fin such 
as might be used in a central fin-tube 
radiator (fig. 1(a)). Only one quar- 
ter of the fin is analyzed (fig. 1(b)) 
because symmetry is assumed about 
x = L and y = 0. There is no heat 
flow across these boundaries. The 


yj Stainless steely t temperature at x = 0 is assumed con- 

1 ss stant and equal to the tube outer- 

surface (base) temperature in a fin- 
tube configuration. 

The analysis is based on steady- 
state two-dimensional conduction in 

(WAS considered in analysis (thickness exaggerated). the fin and the clad with P^fect ther- 

mal contact at the interface (y = t cu J. 

Figure 1. -Cross section of stainless-steel-clad copper fin. The surface of the clad ( y = t ) is 

assumed to be radiating hemispherically 
at constant emittance to space, but radiant interchange with the tubes and in- 
cident radiation from the surroundings are neglected for simplicity. Also, the 
effective sink temperature of space is assumed to be zero. It is assumed that 
the thermal conductivity along the fin does not vary with temperature. For the 
range of variables considered herein the maximum error in assuming constant con- 
ductivity is estimated as less than 0.6 percent by the method of reference 7. 



Basic Equations 


The equation governing the two-dimensional conduction in the copper is 

d 2 T cu U,y) d 2i r cu( x > y) 


and in the stainless steel clad is 

d 2 T ss (x,y) + d 2 T ss (x,y) _ 
dx 2 Sy 2 

where T cu (x,y) is the local temperature in the fin and T ss (x,y) is the local 
temperature in the clad. (All symbols are defined in appendix A.) 

At the boundary y = t cu between the fin and clad, perfect thermal con- 
tact is assumed: 


TcuU^cu) “ ^ss( x ^cu) 
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Also at this boundary there must be continuity of heat flow; hence. 



where k cu is the thermal conductivity of the fin and k ss that of the clad. 
At the radiating surface (y = t), 


-k 


fST ss (x,y)' 


ss 


oy 


Jy=t 


= eoTf s (x,t) 


(5) 


In accord with the assumptions mentioned previously, the equations for the 
remaining boundary conditions are 


T cu U,y) = T ss (x,y) = T 0 
c>T cu (x,y) ^1 ss (x,y) 


dx 


dx 


= 0 


and 


dT cu (x,y) 


w 


= o 


x = 0 


x = L 


> symmetry 


y = 0 


In order to facilitate comparison of the solutions for different cases, 
the following dimensionless parameters are formed: 

(1) Clad-to-core thickness ratio, t ss /t cu 

(2) Aspect ratio of copper fin, t cu /L 

(3) Conductivity ratio, k ss /k cu 

(i) Temperature ratio, 0(X, Y) = T(x,y)/T 0 

(5) Distance ratio for any distance, X = x/L and Y = y/L 

(6) Conductance parameter, A = eaT^l£/k cu t 

This last ratio does not occur naturally in the two-dimensional analysis of 
clad fins, but it is used because of its wide occurrence in fin-tube radiator 
literature . 

Putting the dimensionless parameters into equations (l) and (2) yields 
the following: 
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In the fin: 


In the clad: 


a 2 e cu (x,Y) a 2 e cu (x,Y) _ 
ax 2 + dY 2 

S 2 e ss (X,Y) a 2 0 ss (X,Y) 


ax 2 


dY 2 


= 0 


The boundary conditions then become the following: 
^cu . 


At Y = 


and 


e 4 x ’ ¥) - ¥) 


S0 C1J (X,Y) k ss a0 ss (X,Y) 


= ^ss ss 
k cu 


SY 


W 


At Y = i ^ 

Jj J_l 


1 + 


J ss 


cu 


a0 q .(X,Y) 


_ -y ^cu k cu n 4 / v t\ 

3y A T" k^T 9 ^\ x ’ l ) 


( 6 ) 


(7) 


( 8 ) 


(9) 


( 10 ) 


At X = 0: 


0 ss (X,Y) = 0 cu (X,Y) = 1.0 


( 11 ) 


At X = 1.0: 

a 0 ss (x,Y) a0 cu (x,Y) 
ax ax 


( 12 ) 


At Y = 0: 


a© cu (x,Y) 

aY 


(13) 


It is thus seen that the independent parameters are tq^/L, t ss /t cu , k ss /kcu> 
and A, and their specification determines the temperature distribution through- 
out the fin. Results are given in dimensionless form in terms of these dimen- 
sionless variables; however^ the calculations were made on the dimensional equa- 
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tions (l) and (2). These equations subject to the dimensional boundary condi- 
tions given previously were solved numerically on an IBM 7094 computer using a 
finite-difference, line-by-line, overrelaxation method. Details of the numer- 
ical method of solution are given in appendix B. 


Fin Efficiency 

A common measure of fin performance is the radiating efficiency, which is 
defined as the ratio of the actual heat radiated by the fin based on the calcu- 
lated temperature distribution along the outer surface to the heat radiated by 
an isothermal fin at the base temperature T 0 : that is. 


f ®! s ( x - e ) “ < 14 > 

* / 0 ' 

This parameter will be used to evaluate the effect of various thicknesses of 
clad. 


T] = 


Q 

^ideal 


f ecTf s (x,t) dx 




Effective Thermal Conductivity 

The effective conductivity of a clad fin is a function of the conductiv- 
ities and the relative dimensions of the copper core and the stainless steel 
clad. It is also a function of the path the flow of heat must take to reach 
the radiating surface of the clad fin. These considerations are reflected in 
an effective conductivity k of clad fin which, in this report, is that value 
of conductivity that satisfies 


ecT^L 2 

A eff = — (15) 

where A e ff is that value of A which, when used in one -dimensional radiating 
fin calculations (such as those of refs. 8 and 9 or the constant property cases 
of ref. 7), will result in the value of radiating efficiency T) obtained 
herein by two-dimensional calculations. 


RESULTS AMD DISCUSSION 

Several fin geometries representative of 300-kilowatt and 1-megawatt 
electrical-power-output Rankine-cycle space powerplants were analyzed. Two 
temperature levels (811° and 950° K) were chosen, and the fin length and fin 
thickness were varied to obtain a range of practical values of the identifying 
parameter A. These input parameters are summarized in table I. Each basic 
configuration of table I was analyzed over a range of thickness ratios between 
0 and 1. The variation in thickness ratio for a given configuration was ob- 
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Temperature ratio, TfT ( 


TABLE I. - INPUT PARAMETERS 1 


[Emittance, e, 0.9; thermal conductivity of copper, k c 
3.65 W/(cm)(°K); thermal conductivity of stainless 
steel, k ss , 0.1818 W/(cm)(°K); conductivity ratio, 

k ss/ k ou^ °- 05 -3 


Case 

Base 

Fin 

Copper 

Conductance 

Aspect 


temper- 

length, 

fin 

parameter. 

ratio of 


ature, 

L, 

thickness. 

A 

copper 


V 

°K 

cm 

w 

cm 


fin, 

t /l 

cu' 

A 

811 

2.0 

0.05 

0.1218 

0.0250 

B 



5.1 

.20 

.1904 

.0400 

C 



4.1 

.05 

.4874 

.0125 

D 



10.2 

.20 

.7616 

.0200 

E 



5.6 

.05 

1.0000 

.0089 

F 



15.3 

.20 

1.7136 

.0131 

G 

950 

1.8 

0.05 

0.1499 

0.0286 

H 



4.3 

.20 

.2210 

.0471 

I 



3.6 

.05 

.5996 

.0143 

J 



8.6 

.20 

.8841 

.0236 

K 



12.9 

.20 

1.9892 

.0155 


a To convert from cm to in., multiply by 0.3937; °K to °F, 

°F = 9/5 °K - 460; °K to °R, multiply by 9/5; W/(cm)(°K) to 
Btu/(hr)(ft)(°R), multiply by 57.8. 



tained by decreasing t cu , the 
thickness of copper, and in- 
creasing t ss , the thickness of 
stainless steel clad, while main- 
taining a constant total thick- 
ness t. The conductivity ratio 
k ss /k C u should be a function of 
fin-base temperature; however, 
to reduce the number of parame- 
ters, the same value was used for 
both temperature levels. 

The results to be presented 
are discussed under two main 
headings: first, the degree of 

two dimensionality present in 
the cases analyzed, and second, 
the effect on certain fin parame- 
ters of varying the stainless - 
steel to copper thickness ratio. 


Figure 2. - Surface temperature distribution of cases A and K for thickness 
ratios from 0 to 1.0. _ „ _ 

Degree of Two-Dimensionality 


The items to be discussed 
in this section are temperature 
profiles along the fin length. 
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Distance ratio along fin, X 


temperature drop across the fin 
thickness, effective conductiv- 
ity, and comparison of two- 
dimensional and one-dimensional 
results . 

Temperature drop across fin. - 
Of the several cases analyzed 
(table I), case A has the smallest 
temperature gradient along the fin 
length and case K has the largest. 
The temperature profiles of these 
two extreme cases at three values 
of thickness ratio t ss /t cu are 
shown in figure 2 . The two- 
dimensional effect, that is, the 
difference between the tempera- 
ture profile along the radiating 
surface of the fin (Y = t/L) and 
that along the adiabatic surface 
(Y = 0), is not large enough to 
be seen on the scale of f igure 2 . 


Figure 3. - Temperature ratio drop across fin thickness against dis- The two- dimensional effect 

tance ratio along fin. is shown more specifically on 

figure 3 in the form of the temperature ratio drop across the fin from Y = 0 
to Y = t/L for three values of t ss /t cu . < ^ ie cases shown are case K, because 
it has the steepest X-direction temperature gradient (fig. 2), and case H, be- 
cause it has the largest Y-direction temperature gradient. The variation in 
t ss /t cu vas obtained by reducing t cu and increasing t ss to maintain a con- 
stant total thickness t. The clad fins have a greater two-dimensional effect 
than the all-copper fins (t ss /t cu = 0) as shown in the figure. It should be 
noted, however, that the scale of the ordinate in figure 3 is greatly expanded 
and the maximum difference in temperature ratio between the fin centerline and 
the fin surface is only about 1 . 15 X 10 “3 unit. 


Effective conductivi ty. - As was stated in the INTRODUCTION the value of 
effective conductivity of a radiating clad fin should lie somewhere between the 
one -dimensional value for conduction in the direction parallel to the fin sur- 
face and that for conduction normal to the surface. In view of the small two- 
dimensional effect noted in the previous section, it is to be expected that the 
effective conductivity of the radiating fin would be close to that for one- 
dimensional conduction in the parallel or x-direction. To check this expecta- 
tion, the exact effective conductivity for the two-dimensional radiating fin 
will be investigated and will then be compared with the one-dimensional 
x-direction conductivity. 

An exact value of effective conductivity can be obtained from equation (15). 
Solving_equation (15) for the exact effective conductivity k and forming the 
ratio k/k cu yield 
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Figure 4 - Exact effective conductivity ratio against thick- 
ness ratio. 


cu 


eoT^L 2 

tk cu A eff 


( 16 ) 


The grouping eoToL 2 /tk cu of equation 
(16) is equal to the conductance parame- 
ter A. Substituting A into equa- 
tion (16) yields 


tu 


^eff 


(17) 


A value of A e ff was obtained from 
one -dimensional calculations for each 
value of r\ obtained from the two- 
dimensional calculations reported herein, 
and the effective conductivity ratio was 
obtained from equation (17). On fig- 
ure 4 are shown curves for the exact ra- 
diating-fin effective conductivity ratio 
t ss /t cu for cases E and H. The 

two curves and are omitted 


k/k cu as a function of thickness ratio 
remaining cases listed in table I fall between these 
for clarity. The close agreement among the several cases indicates that the 
effective conductivity ratio k/k cu is dependent primarily on thickness ratio 
ts S /tcu for a given conductivity ratio k ss /k cu . The exact value of effective 
conductivity will be compared with the effective conductivity based on one- 
dimensional conduction in the x-direction, which is the thickness-averaged 
value given by 




^SS^-SS + ^CU^-CU 

^ss + ^cu 


(18) 


Writing equation (18) in terms of dimensionless ratios yields 

_ ^ss k ss + ^ 

k x _ ^cu k cu 

]r ~~ +• 

CU ^ + °SS 

^CU 


(19) 


Thus, the ratio k x /k cu is a function only of the thickness ratio 1 SS A CU 
and the conductivity ratio k ss /k cu . The relation between k x /k cu and tss/^cu 
is y shown on figure 5(a) for a conductivity ratio k ss /k cu of 0.05. Also 
plotted in figure 5(a) is the exact effective conductivity ratio of the case 
that deviates most from the thickness -averaged conductivity (case H). The close 
agreement between k x /k cu and k/k cu is evident. 


As was stated earlier, only one value of k ss /k GU was used even though 
two temperature levels were considered. To determine whether the change in 
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(a) Conductivity ratio, 0.05. (b) Conductivity ratio, 0.07. 

Figure 5. - Comparison of thickness-averaged conductivity with exact two-dimensional effective conductivity. 


k ss /k cu over the temperature range of interest significantly affected the 
agreement of k/k cu and k x /k C uj case H was analyzed with a conductivity ratio 
k ss /k cu of 0,07 (corresponding to a temperature of 911° K). The results are 
shown in figure 5(b), and the agreement is seen to be not significantly differ- 
ent from that of figure 5(a), 

This close agreement of k x /k cu with k/k cu means that the thickness- 
averaged effective conductivity can safely be used in radiator heat-transfer 
calculations for geometries comparable to those treated herein. 

Comparison of one -dimensional and t wo-dimensional results . - As a further 
check of the closeness between the results of one-dimensional calculations with 
the thickness-averaged conductivity k x and the two-dimensional solutions, 
several cases were analyzed by the method of reference 7 but with constant 
properties. (The one -dimensional methods of refs. 8 or 9 could also be used.) 
The cases chosen from table I were those having extreme values of A, cases A 
and K, and those having extreme values of t/L, cases E and H (case H has the 
greatest two-dimensional effect). In all cases the temperature profiles and 
fin efficiencies agreed with the two-dimensional results to within three places 
with small variations in the fourth place. Thus, for the range of A, t/L, and 
t ss /t cu treated herein, the use of thickness -averaged conductivity in one- 
dimensional calculations will give essentially the same results as the more 
exact two-dimensional calculations. 


Effect of Varying Thickness Ratio 

The stainless-steel to copper thickness ratio used will in general depend 
on the particular application. All the cases of table I were analyzed with 
thickness ratios t ss /t cu of 0, 0.1, 0.3, 0.5, and 1.0. The effect of this 
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Heat rejection per unit weight ratio, (Q/W)/(Q/W) ( 



Conductance 
_ parameter, . 


Clad-to-core thickness ratio, t ss /t cu 

- Fin efficiency against thickness ratio for cases with constant 
i thickness. 



Conductance 
. parameter, 


.2 .4 .6 .8 1.0 

Clad-to-core thickness ratio, t ss /t cu 

Figure 7. - Heat rejection per unit weight ratio against thickness ratio. 


variation m thicknsss ratio on 
several parameters used in evalu- 
ating fins or fin materials is 
presented in this section. Most 
of the results presented were 
obtained from the same two- 
dimensional calculations as the 
temperature profiles presented 
previously. Essentially the 
same results could also have 
been obtained from one- 
dimensional calculations because 
of the close agreement noted in 
the preceding section. 

Fin efficiency . - The effect 
on fin efficiency tj of varying 
^ss/^cu shown on figure 6. 
This figure indicates a signifi- 
cant drop in heat transfer as 
the high conductivity copper is 
replaced by the lower conductiv- 
ity steel. Thus., making the 
simplest assumption and treating 
the clad as if it were copper 
would introduce a large heat re- 
jection error in radiator calcu- 
lations (except for very small 
tss/tcu)- Fo r example, for case K 
(A = 1.99) with t ss /t cu = 0.5, 
the fin efficiency is 16 percent 
lower than that for the all 
copper fin (t ss /t cu =0). How- 
ever, as was indicated earlier, 
the clad fin can be treated 
accurately in one -dimensional 
radiator calculations with the 
use of the thickness -averaged 
( effective ) conductivity. 

H eat rejection rate per unit 
weight . - An important parameter 
used for comparative evaluation 
of a radiator component is the 
heat rejection rate per unit 
weight Q/W. Therefore, to study 
the effect of varying the propor- 
tion of stainless steel to copper 
in a clad fin the ratio 
( Q/w )/ (Q/w) is formed, where 


for an all 


(Q/W) 0 is the Q/W for an all 
copper fin (t ss /t cu = 0). This 




Figure 8. - Ratio of effective conductivity to effective density 
against thickness ratio. 


stainless-steel-clad copper on this bai 
is the average density of the clad fin 


ratio is plotted in figure 7 as a 
function of t SB /t cu . For most practi- 
cal radiators , t ss /t cu will be less 
than 0.5, and consequently the reduc- 
tion in Q/W will vary between 3 
and 14 percent for values of A 
ranging between 0.12 and 2.0, respec- 
tively. However, this is the effect 
on the fin only. The net effect of 
using clad fins on the total radiator 
Q/W will depend on the relative con- 
tributions of the fins to total radia- 
tor weight and in many cases will be 
less than the effect shown in figure 7. 

Conductivity to density ratio . - A 
parameter based on material properties 
that is used to compare the conducting 
capabilities of materials on a weight 
basis is k/p_(ref. 3). To compare _ 
is, the ratio k x /p" is formed, where p 
md is given by 


P 


^ss^ss + ^cuPcu 
"^ss + ^cu 


( 20 ) 


and k x is the thickness -averaged conductivity defined previously in equa- 
tion ( 18 ) . 


TABUS II. - CONBUCTIV ITT TO DENSITY RATIO a FOR 
SEVERAL MATERIALS AT 811° K 


Material 

Conductivity to 
density ratio, 

k x / p 

(W)(m 2 )/(kg)(°K) 

Pyrolytic graphite 

0.132 

Beryllium 

.0583 

Commercial graphite 

.0550 

Copper, t ss /t cu = 0 

.0428 

Stainless -steel-clad copper, t ss /t cu = 0.5 

.0308 

Stainless-steel-clad copper, t ss /t cu = 1.0 

.0231 

Molybdenum 

.0132 

Columbium 

.0077 

Tantalum 

. 00385 

Titanium 

.00385 

Stainless steel 

.00286 

Vanadium 

.00231 


a All data except stainless-steel-clad copper obtained from 
fig. 9 of ref. 3. 
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The ratio S^/p is plotted against t ss /t cu in figure 8. In table II 

these values of k x /p" are compared with the k/p values obtained for a variety 
of materials from reference 3. It can be seen that for thickness ratios up to 

0.5, stainless-steel-clad copper has a higher k/p" value than most of the 
materials listed, and is not too far removed from the best materials . 


SUMMARY OF RESULTS 

For the particular configurations analyzed the following results have been 
obtained: 

1. The temperature in the clad fin shows only a slight two-dimensional 
effect. 


2. The effective conductivity of a radiating clad fin is very nearly equal 
to the thickness-averaged conductivity (one-dimensional) of the clad fin. 

3. One -dimensional calculations that use the thickness -averaged conductiv- 
ity yield essentially the same results for temperature profile and fin effec- 
tiveness as the two-dimensional calculations. 

4. Increasing the proportion of stainless steel in a clad fin of given 
total thickness decreases the radiating effectiveness of the fin. For a 
stainless-steel to copper thickness ratio of 0.5, the reduction is 16 percent. 

5. The ratio of effective conductivity to effective density of a stainless- 
steel-clad copper fin is relatively high and remains competitive with the other 
materials of interest. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, August 17, 1965. 
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APPENDIX A 


s4 

a 

86 

b 


c 

a 

i 

j 

k 

k 

L 

m 

NDX 

NDYC 

NDYF 

NX 

NY 

n 

Q 

^ideal 


SYMBOLS 

matrix of system of finite-difference equations 
diagonal elements of matrix 86 

submatrices of order NY of d which lie along the diagonal 
off-diagonal elements of matrix 
off-diagonal submatrix of order Wf of s/ 
elements of matrix # 

vector whose components are right-hand sides of finite-difference 
equation 

position number of horizontal grid line 
position number of vertical grid line 
thermal conductivity, W/(cm)(°K) 

effective conductivity of radiating clad fin, w/(cm)(°K) 

effective conductivity of clad fin transferring heat by conduction in 
x-direction only, W/(cm)(°K) 

fin length, cm 

number of current iteration 

number of vertical grid strips 

number of horizontal grid strips in clad 

number of horizontal grid strips in fin 

total number of vertical grid lines 

total number of horizontal grid lines 

distance normal to surface, cm 

heat rejection rate, W/(cm)(sec) 

heat rejection rate of an isothermal fin at temperature T Q , 
w/(cm)(sec) 
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t EH *TEH 


s distance along a surface, cm 

T temperature, °K 

vector of temperatures 

vector of temperatures, intermediate in the computation 
t thickness, cm 

W fin weight, kg 

X x/L 

x distance along length of fin, cm 

Ax grid spacing in x-direction, cm 

Y y/L 

y distance along thickness of fin, cm 

grid spacing in y-direction, cm 
"z dummy variable vector 

components of z 

coupling parameters (eg.. (B15)) in difference equations 
e emittance 

T] radiating efficiency of fin 

6 temperature ratio, T/To 

3 2/ 

A conductance parameter, €gT 0 L /k cu t 

A e ff A of clad fin, eaT^L 2 /kt 

p vector of ratio of elements of the matrix & 

p^ components of p 

p dens i ty , kg/ m 3 

p* thickness -averaged density, kg/m 3 

o Stefan-Boltzmann constant, 5.67X10“^ W/(m 2 )( 0 K^) 

cd block-successive overrelaxation parameter 



optimum value of cd 


®b 

Subscripts : 
cu copper fin 

J variable horizontal position along fin 

NX last horizontal position along fin 

NY last vertical position on fin 

o value at x = 0 or t ss /t cu = 0 

ss stainless steel clad 

Superscripts : 

(m) m^* 1 iterate 

(o) initial approximation 
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APPENDIX B 


NUMERICAL APPROXIMATION OF TEMPERATURE DISTRIBUTION 
Equations (l) and (2) my be replaced by the single equation 


(Bl) 


vhich holds in the entire fin as a partial differential equation with discon- 
tinuous coefficients. Such problems have been treated recently (e.g., ref. 10). 
To solve equation (Bl) subject to the boundary conditions , a rectangular grid 
is superimposed on the composite rectangle above the x-axis of figure l(b) (p. 3) 
in such a way that each boundary (including the interface) is a grid line, Ax 
is a constant L/NDX where NDX is the number of vertical strips laid off in 
the rectangle by the discretization, and Ay is one constant in the fin 
Ay cu = tcu/NDYF and another in the clad Ay ss = t ss /NDYC. The grid points are 
labeled with ordered pairs (l,J) of positive integers in the manner to be de- 
scribed. 


Let NY = NDYF + NDYC +1 be the total number of horizontals in the grid. 
Number these 1, 2 ... NT starting at the lowermost (y=0). Let NX = NDX + 1 
be the number of verticals in the grid. Number these 1, 2 ... NX from left 
to right starting at the leftmost (x = 0). The first element of the ordered 
pair (IjJ) is the position number of the horizontal and the second element is 
that of the vertical; thus, (l,J) is such that 1 < I < NY and 1 < J < NX. 

Next the partial differential equations and boundary conditions are re- 
placed by a system of algebraic equations (some of which are nonlinear) ob- 
tained by integrating the partial differential equation over the cells of the 
"dual grid" (ref. 11, p. ll), by applying the divergence theorem, and by re- 
placing the normal derivatives occurring in the boundary integrals by differ- 
ence quotients, zero or -eaT 4 , depending on the position of the cell in ques- 
tion (ref. 12, sec. 6.3). 


This process is illustrated by the following example taken at the point 
x = Ax, y = t (i.e., (l,J) = (NY,2)). The cell in question of the dual grid 
(shaded in fig. 9) has its boundaries (except the upper one) midway between 
grid lines of the original grid. The algebraic equation at (NY, 2) is obtained 
in the following steps. 


0 = 


ff [s ( k si) + ir 

Cell 



dx dy 


0 = 


£ 


Around 

cell 



ds 
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Figure 9. - Detail of grid for numerical solution. 


0 - k, 


SS 


~T 0 - T(NY, 2 )" 

A^ss + k 

T(NY - 1,2) - T(NY,2l 

Ax 

2 ss 

A^ss 


+ k 


ss 


T(UY ; 3) - T(NY,g 2 ffiss _ eor^jjy^Ax (B2) 
Ax 2 


This particular point was chosen to illustrate the fact that sometimes a path 
length in the integration is only one-half a total mesh width (hy ss /2) and to 
illustrate the ease of handling the two types of boundary conditions inherent 
in this method of discretization (note first and fourth terms). 

Repeating the integration on all cells that surround the points ( I, J) , 

1 < I < NY and 2 < J < NX, yields a system of NY(NX - l) equations (of which 
NX - 1 are nonlinear, viz., (NY,J), 2 < J < NX) in as many unknowns T(I,J). 

It is noted that T(l,l) = T 0 when 1 < I <1 NY. 

A direct solution of this system is impossible because of the nonlinearity; 
therefore, an iterative procedure must be used. Let T(m) (I,J) represent the 
m th iterate for T(I,J) where T(°) (I,j) is some initial approximation to 
T( I, j) . Equations at the points (NY,J), 2 < J < NX, are linearized in the same 
way that Newton's method for finding the roots of equations is derived (see 
ref. 13, p. 192). The term [T( m )(NY,J)] 4 is replaced by the first two terms 
in its Taylor expansion about T( m- ^- )(NY, J), that is by 

[ T (m-1)(NY,J)] 4 + 4[T( m - 1 )(NY,J)] 3 [T( m )(NY,J) - T( m ~ 1 )(NY, J)] 


which reduces to 


4[T(m-l)(NY,J)] 3 T (m)(NY,J) - 3 [t( 111 - 1 )(NY, J)] 4 

The resulting system of equations linear in T( m ) can be written in matrix nota- 
tion as 

sf T^) = d (B3) 


where 

y = ^[5(m-D] (B4) 

is a square matrix of order NY(NX - l) and 
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1 - a(5(«-D) 


( B5 ) 


is a column vector with NY(NX - l) components. is a block tridiagonal 

matrix, that is. 




Here 


[*z 

-V 

0 . 

• . 0 



^3 

-<e . 

. . 0 

0 \ 

0 

• 

-Sf 

• 

#4 . 

• • 

. . 0 

° 

• 

• 

• 

: 0 

• 

• 

0 

• • 

• • 

0 . 

• » • 

* • • 

• * ^NX-1 

-v I 

v° 

0 

0 . 


&WX J 

^( T (m- 

■ 1 )( 

NY,J)) 

2 < J < 

NX - 


(B6) 


(B7) 


are square, tridiagonal matrices of order ^ NY, all of whose elements are con- 
stant (independent of both J and T( m- l)) except for the element in the NY, NY 
position which depends on (NY, J) . The remaining matrix on the diagonal 

of sd is 


#NX = | (NY, NX)) (B8) 

and ^ is a diagonal matrix (constant) of order NY. The iterative method may 
be described as a line-by-line (vertical) overrelaxation scheme, which is 
carried out in the following way. 


An initial approximation to the temperature is taken at all points (e.g., 
t(°)(I,J) = T Q ). The temperatures are changed by sweeping through them from 
left to right an entire vertical at one time; the sweep is repeated if the con- 
vergence criterion has not been met. If attention is focused on the J ^ verti- 
cal, the change is carried out by the following formulas . Let 



i(m) 


(1,J) 


T (m)(NY,j) 


(B9) 


be the vector of temperatures along the vertical at the iterate. It is 
assumed that T^- m ^ has already been computed and T^. m ^ defined by 



( Bio) 
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Start 



Figure 10. - Computer flow diagram. 


vhere % is stored in computer 
memory in partially factored form 
since is tridiagonal and 

almost constant.^ Therefore the 
computation of Tj is only a mat- 
ter of a correction to the right- 
hand side and a back substitution. 
Then, 

7p(m) _ ^(m-l) + ^r^m) _ 
j J L J J J 

(BID 

where cd is the "block successive 
overrelaxation parameter" (ref. 12, 
ch. 4) that controls the rate of 
convergence, and 

0 < cd < 2 (B12) 

To further clarify this, the 
computation (eq.. (BIO)) on the 
second vertical (J = 2) is examined 
in detail, and reference is made 
to the block diagram of the com- 
puter program ( fig . 10 ) . 


For J = 2, 


and 


'2 = 


a l 

1 

cr 1 

H 

0 ... 

0 

0 

0 

u 1 

H 

a 2 

" b 2 • ’ * 

0 

0 

0 

0 

“ b 2 

a 3 ..• 

0 

0 

0 

. 

. 

• • • • 

• 

• 

• 

• 

• 

• • • • 

• 

• 

• 

• 

• 

• • • • 

• 

• 

• 

0 

0 

0 ... 

a WY-2 

-b EY-2 

0 

0 

0 

0 ... 

_b NY-2 

a NY-l 

-b NY-l 

0 

0 

0 ... 

0 

" b NY-l 



*€ = diagCcj . . . Cjjy ) 


(B13) 


(B14) 
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For convenience, define 


an — 


^ss^ss 


Ax 


cu^cu 
a 2 Ax 


Axk 


ss 


°° 3 A^ss 


Axk, 


CL A 


cu 




cu 




(B15) 


Then, when equation (B2) and analogous equations for the other points on 
this vertical are considered, the elements of &2 an( l ^ can Be written as 
follows : 


< 


a 2 + 
2(a 2 
+ 

2(aQ_ 

k. a l + 


a 4 

+ a 4 ) 

a 2 + a 3 + a 4 
+ a 3 ) 

a 3 + 4ec Ax[T( m - 1 )(WY,J)] 3 


i = 1 

i « 2 . . . NDYF 
i = NDYF + 1 

i = NDYF + 2 . . . NY - 1 
i = NY 


b. 


1 


fa 4 

\a 3 


i = 1 . . . NDYF - 1 
i = NDYF ... NY - 1 


r ( B16 ) 


c 


i 


< 


a 2 /2 

a 2 

( an + a 2 )/z 


a l 

^0-1/2 


i = 1 

i = 2 . . . NDYF 
i = NDYF + 1 

i = NDYF + 2 . . . NY - 1 
i = NY 


Therefore, with equation (BIO) rewritten for J = 2 as 


# 2 % m ) =<r[T^ m ) + 

the computation for T^ m ^ is carried out by solving the following system of 
equations where, for brevity, the vector 
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(B17) 



replaces the vector T 2 m ^ : 


a l z l “ b 2 z 2 
~ b l z l + a 2 z 2 ~ h 2 z 3 
- b 2 z 2 + a 3 z 3 


h 3 z 4 


Cl [T 0 + T(l,3)] 
c 2 [T q + T(2,3)] 
c 3 [T q + T(3,3)] 


y (bib) 


-' b NY-2 z NY-2 + a M-i z NY-i “ b NY-i z mr 
- %£-l z M-l + a NY z WY 


c m-it T o + t(my-i,3)] 
c m [T 0 + T(nr,3)] 


Here the superscripts on T(l,J) have been dropped to be consistent with what 
actually takes place in the computer program, namely, the immediate replacement 

of tO 31 " 1 ) by once it is computed. 


Of all the numbers aq, bq and cq, only aj^y depends on the fact that 
computations are being carried out on the second vertical and at the m zn 
iterate. Thus, if some other vertical, other than the last, were examined, the 
only changes that would occur in the coefficients in (B18) are as follows: 

(1) ajjy is replaced by its correct value. 

(2) T 0 + T(I,3) is replaced by T(l,J - l) + T(I,J + l), I = 1 . . . NY. 

Therefore, if the previous system is solved by the elimination of unknowns, it 
is seen that as one proceeds from vertical to vertical and from iteration to 
iteration, the same arithmetic operations are performed on the same members with 
the same ratios (pi) appearing repeatedly (at least so far as the left-hand 
sides are concerned). Even for the right-hand sides the same operations are 
being performed with the same ratios (pq) as with the left-hand sides, but the 
operations are performed on numbers that depend on J and m. It is therefore 
advantageous to compute and store these pq ratios before any iterations are 
performed. It is this saving of the p^ along with the modified aq that has 
been referred to previously by "triangularizing of $g<^' (block 2 of fig. 10). 
Applying the previously mentioned operations and p^ ratios to the iteration- 
and vertical-dependent right-hand sides (blocks 4 and 5 of fig. 10 ) leads to a 
reduced system of equations, the last of which involves only zj^y (block 6). 

This is then substituted into the preceding (modified) equation, which involves 
it with ZKjy^p and the process continues backwards (in I) to zq (blocks 7 
and 8). The values for the temperatures on this second vertical are then re- 
placed by the new values (block 9): 
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T(I,2) + cd [ zj - T(I,2)] 


1=1. . .NY 


Numerical experiments at Lewis indicate that there is an optimum value 
for ay, as in the case of equations that are all linear , where an a priori 
value can be given. For example , for a case where NDYF = 4, NDYC = 2, 

NDX = 100, L = 0.1333, t cu = t ss = 0.0004167, k ss = k cu = 210, T 0 = 2337, 
e = 0.9, and a = 0.1713X10“®, it is found that cd^ near 1.96 gives a relative 
residual error 

| heat in - heat out I 
heat out 

near 0.0001 in less than 300 iterations. This is a practicable convergence 
criterion for the problem at hand and is so used in the computer program 
(block 11 of fig. 10). 

Here 

k ^ ds 
on 

end 


heat in = 


-JL 


Left 


heat out = 


L 


Radiating 

boundary 



ds 



ds 


Trapezoidal integration is used to be consistent with the finite-difference 
equations. These experiments also indicate that for cd < heat in increases 
monotonically to meet heat out, which is decreasing as m increases if 

is chosen to be T Q everywhere. 

The program was converted to double precision after it was found that the 
large aspect ratio of the cells of the grid led to appreciable loss of signifi- 
cance. The program takes less than 1 minute to carry out 400 iterations for 
cases similar to those mentioned previously. 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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